Transmission of matter wave solitons through nonlinear traps and barriers 



O 
O 

(N 

X 



Josselin GarnierfTT and Fatkhulla Kh. AbdullaevJ 

f Laboratoire de Probabilites et Modeles Aleatoires & Laboratoire Jacques-Louis Lions, 
Unwersite Paris VII, 2 Place Jussieu, 75251 Paris Cedex 5, France 
J Physical- Technical Institute of the Uzbekistan Academy of Sciences, 
700084, Tashkent-84, G.Mavlyanov str.,2-b, Uzbekistan 

The transmissions of matter wave solitons through linear and nonlinear inhomogeneities induced 
by the spatial variations of the trap and the scattering length in Bose-Einstein condensates are 
investigated. New phenomena, such as the enhanced transmission of a soliton through a linear trap 
by a modulation of the scattering length, are exhibited. The theory is based on the perturbed Inverse 
Scattering Transform for solitons, and we show that radiation effects are important. Numerical 
simulations of the Gross-Pitaevskii equation confirm the theoretical predictions. 
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I. INTRODUCTION 

The transmission of matter wavepackets through inho- 
mogeneities of different types of Bose-Einstein conden- 
sates (BECs) has recently attracted a lot of attention, be- 
cause this phenomenon is important for the design of con- 
trol methods of the soliton parameters and atomic soliton 
lasers 0] . The transmission and reflection of bright and 
dark matter wave solitons has been studied in the case 
of linear inhomogeneities, induced by the variations in 
space of the potential field 0, H S @. The effect of 
a potential step or impurity, including the soliton train 
evolution, has been analyzed in Q . The case of inhomo- 
geneities produced by spatial variations of the scattering 
length has been less investigated 0, @, 0- Numerical 
simulations of the problem of soliton propagation have 
been carried out, when the linear and nonlinear inhomo- 
geneities compete with each other. It was found that for 
some choices of parameters the enhanced transmission of 
soliton through the nonlinear barrier is possible |lfj| . 

The purpose of this work is to develop the theory de- 
scribing the transmission of matter wave solitons through 
nonlinear barriers. Such barriers can be produced by us- 
ing the Feshbach resonance method, namely by the lo- 
cal variation of the external magnetic field B(x) in space 
near the resonant value B c 0] . By the small variation of 
the field near the resonant value we can induce the large 
variations of the scattering length in space according to 
the formula 



a s (x) = ttb I 1 - 



B c - B{x) 



where a& is the background value of the atomic scattering 
length and A is the resonance width. Optical methods for 
manipulating the value of the scattering length are also 
possible As a result, the mean field nonlinear coeffi- 
cient (which is proportional to the scattering length a s ) in 



the Gross-Pitaevskii equation has a spatial dependence. 
We will use the perturbed Inverse Scattering Transform 
theory (see for example 003111) 

to describe the trans- 
mission of bright matter wave solitons through the non- 
linear barriers. This approach allows us to analyze the 
adiabatic dynamics of solitons as well as the radiative 
processes during the soliton propagation through inho- 
mogeneities. 



II. THE MODEL 

The quasi one-dimensional Gross-Pitaevskii (GP) 
equation describing the wavefunction of BEC in an elon- 
gated trap has the form [Tj| 



h 2 

ih4>t + -z—^xx 
2m 



V(x)iP-g 1D (x)\iP\ 2 iP = 0. (1) 



Here gio(x) — 2hu>j_a s (x), where lo±_ is the transverse os- 
cillator frequency, a s (x) is the spatially dependent atomic 
scattering length, and V(x) is the linear potential. Both 
a s {x) and V{x) are assumed to be constant outside 
a given domain, where we assume that the scattering 
length takes the constant negative value a s o and V is 
zero, J dx = N, where N is the number of atoms. 
We denote by g\o = 2Huj±a s o the reference value of the 
nonlinear coefficient. 

From now on we express x in units of the healing length 
£ = h/y/nogiDm and t in units of t,Q = £/(2c), where 
c = nogm/m is the Bogoliubov speed of sound and no 
the peak density. Wc also normalize the mean field wave- 
function by ^/no, so that we obtain the dimcnsionlcss GP 
equation in the form of the Nonlinear Schrodinger (NLS) 
equation 



iu t + u xx + j{x)\u\ 2 u - Vi(x)u = , 



(2) 



with j(x) = 2 - V nt (x) = 2a s (x)/a s0 . If Vi = V n i = 
0, then this equation can be reduced to the standard 
NLS equation that supports soliton solutions. The bright 
soliton solution is: 
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u s (x,t) = 2z/sech[2j/(x— x s (t))] exTp[2i/i(x—x s (t))+i(j) s (t)]. 
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The soliton amplitude is 2v and its velocity is 4/x. The 
soliton center and phase x s (t) and <f> s it) satisfy 

The matter wave soliton moving in the linear and nonlin- 
ear potentials Vi and V n i experiences velocity and mass 
modulations and emits radiation. To describe this pro- 
cess we use the perturbation theory based on the Inverse 
Scattering Transform (1ST). 

III. QUASI-PARTICLE APPROACH 

Applying the first-order perturbed 1ST theory, we ob- 
tain the system of equations for the soliton amplitude 
and velocity 

dv da 1 ,. dx s 

where the prime stands for a derivative with respect to 
x and the effective potential has the form 

W(u,x) = Wi(u,x) + W n i(v,x), 

with 

f°° 1 z 
Wi(y,x)=v V l (—+x)dz, (3) 

J-oo coslr(z) 2v 

W nl {v,x) = 2v 3 \ Vnl (—+x)dz. (4) 

J-oo cosh (z) 

We can thus write the effective equation describing the 
dynamics of the soliton center as a quasi-particle moving 
in the effective potential W: 

v^ = -W'{v Q ,x s ), (5) 

where 4fo is the mass (number of atoms) of the incoming 
soliton. This system has the integral of motion 

¥ (jt) 2 + w{y Q ,x s ) = Sv^l (6) 

where is the velocity of the incoming soliton. Note 
that this approach is standard and it was recently ap- 
plied in In this work, it is the first step as we will 



include second-order and radiation effects in the next sec- 
tion. Let us briefly discuss the main results that can be 
obtained with the quasi-particle approach. 

Barrier potential. Let us first examine the case 
where the potential V is a barrier, meaning that V > 
and lim^i^oo V(x) = 0. When the soliton approaches 
the barrier, it slows down, and it eventually goes through 
the barrier if its input energy is above the maximal en- 
ergy barrier, meaning 8^0/^0 > W max {vo). For a lin- 
ear barrier potential which is an even function, we have 
W max (v) = Wi(f, 0), where Wi is given by |J3J|. For a 
nonlinear barrier potential which is an even function, we 
have W max {v) — W n i{v, 0). After passing through the 
barrier, the soliton recovers its initial mass and velocity. 
In that sense, the transmission coefficient is one. 

If, on the contrary, the velocity of the incoming 
velocity is such that 8fo/io < W max (vo), then the soliton 
is reflected by the barrier. After the interaction with 
the barrier, the soliton velocity takes the value — 4/io- In 
that sense, the transmission coefficient is zero. 

Trap potential. We now examine the case where 
the potential is a trap, meaning that V < and 
limui^oo V(x) — 0. When the soliton approaches the 
trap, it speeds up, and it eventually goes through the 
barrier whatever its initial velocity is. This is the pre- 
diction of the quasi-particle approach. However, we shall 
see that the interaction with the trap generates radiation 
and reduces the mass and energy of the soliton. As a re- 
sult, the soliton may not be able to escape the trap if its 
initial velocity is too small. We shall discuss this point 
in the next section. 



IV. RADIATION EFFECTS 

The properties of the radiation emitted by the soli- 
ton interacting with the inhomogcneities are defined 
by the Jost coefficients a(t, A) and b(t, A) of the as- 
sociated linear spectral problem for the NLS equation 
[T3L [l4l ITEf . We assume that the linear and nonlin- 
ear potentials are localized functions, and we denote by 
V(k) — J V(x) exp(ikx)dx their Fourier transforms. If 
the soliton goes through the potential, then we find 



-f> = +00, A) 
a 



—iir(\ — fi— ivy 



16/i 3 cosh ^ 



7T \ 2 +V 2 -H 2 

4 nu 



[(A + M ) 2 + ^ 2 ][(A-3 M ) 2 + 8 M 2 + 



l -V nl (fc(A (I /,/i)) 



where 



(a-m) 2 + 



The total radiated mass density during the interaction 
with the potentials is 

2 



n(A) 



-(t = +00, A) 
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The total mass (number of atoms) and Hamiltonian are 
preserved: 



=0.5 



TV ■ 



|it| 2 diE 



H = / [\u x \ 2 - \u\ 4 + Vi(x)\u\ 2 + -V nl (x)\u\ 4 ]dx. 

They can be expressed in terms of the radiation and soli- 
ton components as 



N = Av+ I n(X)dX, 

J — oo 

16 f°° 
H = 16^/i 2 - — v 3 + 2W(u, x s )+4 / X 2 n(X)dX. 

•J J -oo 

Note that this expression of the Hamiltonian generalizes 
the expression of the energy ©. The coefficients (vt, Mt) 
of the transmitted soliton are 



1 



n(X)dX, 



(7) 



= A^o - a / + n{X)dX. 8) 

These formulas allow us to study and characterize the 
transmission of a soliton through a general barrier in vari- 
ous regimes. We can consider the transmission/reflection 
of a bright soliton through a nonlinear barrier, the trans- 
mission/trapping in a nonlinear trap, and competition 
effects for the transmission through the superposition of 
nonlinear and linear potentials. 

In the next subsections, we compare our theoretical 
predictions with results from numerical simulations of 
the one-dimensional GP equation, with a Gaussian linear 
potential and/or a nonlinear Gaussian potential: 



Vi (x) = Vio exp - -r 



V n l (x) = V nW exp -■ 



We consider the transmission of a soliton incoming 
from the left homogeneous half-space with parameters 
(vq,IJ,o). We plot the variations of the soliton param- 
eters versus the value of the velocity for vq = 0.5 and 
x c = 0.5. We consider different combinations of linear 
and nonlinear potentials. 

Nonlinear barrier. We consider the case Vq = and 
Vnio > 0. We have seen in Section ITTll that the condition 
for the transmission through a nonlinear barrier is that 
the velocity of the incoming velocity is large enough so 
that 8^oMo > W max (po). Note that this means that the 
critical velocity parameter ^, cr it, defined by the identity 
%v$\i\ rit = W max (p ), is of the order of V^ 2 for V nW < 
1. If the transmission condition fails, then the soliton is 
reflected. 

These results are obtained in the quasi-particle ap- 
proach and neglect the radiation emission phenomemon. 
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FIG. 1: Soliton parameters variations versus input velocity 
parameter fio for a nonlinear barrier. The theoretical reflec- 
tion condition is 8^o/io < Wmaxi^o) which gives \i C rit = 0.21 
for V n io = 1 and /j, C rit = 0.15 for V„io = 0.5, in excellent 
agreement with the simulations. The predicted mass and ve- 
locity reductions are also quantitatively accurate, especially 
for V n io = 0.5. 



Taking into account radiation yields that the transmis- 
sion is not complete in the case 8fo/Xg > W max (vo), in 
the sense that the transmitted soliton mass is not equal 
to the incoming soliton mass (see Figure^. The mass 
loss is described by (0) , and it is of the order of V n io for 
Vnio < 1. 

Nonlinear trap. We consider the case Vio = and 
Vnio < 0. The quasi-particle approach predicts full 
soliton transmission, but neglects radiation phenomena. 
When taking into account radiation emission, we can ex- 
hibit the mass and energy loss during the interaction with 
the trap potential. The losses are described by (|7I8|I and 
are accurate for an initial velocity large enough. Indeed, 
if the initial velocity is not large enough, then the energy 
loss during the interaction with the potential does not 
allow the soliton to escape the trap. 

For small initial velocity, the expression for the radi- 
ation is not precise enough, because the velocity experi- 
ences a modulation which is relatively large compared to 
its initial value. When the soliton center is x, the soliton 
velocity is given by 4/i(a;) with 



H(x) 



W n i{vQ,x) 



16^o 



Its maximal value is [i max = vVo W n i{v , 0)/(16^ ) 
where W n i(vo, x) < is given by During the interac- 
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FIG. 2: Soliton parameters variations versus input velocity 
parameter fio for a nonlinear trap. The theoretical trap con- 
dition is ht < which gives [i C rit ~ 0.12 for V n io = — 1 and 
Merit = 0.07 for V n io = —0.5. The agreement with the simula- 
tions is noticeable for V n io = —0.5, while the case V n io = — 1 
is only in qualitative agreement with the simulations. 



FIG. 3: Soliton parameters variations versus input velocity 
parameter /io for a linear trap. The theoretical trap condition 
is /it < which gives fi C rit = 0.15 for Vio = — 1 and fi cr it = 
0.1 for Vio = —0.5, in good agreement with the simulations in 
the case Vio = —0.5. 



tion with the barrier, the soliton velocity is about 4/Li max , 
so the energy loss can be estimated by —4 J X 2 n(X)dX, 
where [i max is substituted for mo in the expression of 
b/a. If a negative value mt is obtained in @, then this 
means that the soliton has not been transmitted, but 
was trapped by the nonlinear potential (see Figure 
Note that n(A) ~ V^ lQ so that the critical value Mcri* for 
trapping is of the order of V n io- 

If the initial velocity is large enough to ensure trans- 
mission, then the soliton emits radiation and looses 
mass. The mass loss is described by Q and it is of the 
order of V n io for V nm < 1. 

Enhanced transmission by nonlinear modula- 
tion. As pointed out in |17| . the nonlinear potential 
V n i can help a soliton going through a trap potential Vi . 
We illustrate this assertion based on numerical experi- 
ments in this section and justify it with our perturbed 
1ST approach. 

By comparing the transmission through a linear trap in 
presence or in absence of a nonlinear positive potential, 
we confirm the numerical conjecture that the transmis- 
sion coefficient can be significantly increased by a nonlin- 
ear modulation (see Figures |3I4|) . In fact, the radiation 
emitted by the soliton due to the interaction with the lin- 
ear trap and with the nonlinear potential can cancel each 
other, resulting in an enhanced soliton transmittivity. 

A similar result can be obtained with a linear barrier. 
A nonlinear negative modulation can help the soliton go- 
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FIG. 4: Soliton parameters variations versus input velocity 
parameter fio for a linear trap superposed with a nonlinear 
barrier. The theoretical trap condition is [Xt < which gives 
Merit = 0.08 for V n i = 1, V i0 = -1 and Merit = 0.04 for Vnio = 
0.5,V;o = —0.5, in good agreement with the simulations in the 
case Vnio = 0.5, Vio = -0.5,. 
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ing through the linear barrier by reducing the radiation 
emission and by reducing the minimal velocity reached 
by the soliton during the interaction. This was predicted 
by numerical simulations in |10| . 

For consistency we propose an experimental config- 
uration where the enhanced transmission could be ob- 
served for the 7 Li condensate. The transverse frequency 
u>± w 2-7T x 10 3 H and the density is n = 10 9 in -1 . 
The healing length and speed of sound are £ ~ 2 fxxn 
and c ~ 5 mm/s, respectively. As a typical experiment, 
we could consider a soliton with about 10 3 atoms and 
width s» 2£ rj 4 /im (so that v — 0.5) and a linear 
trap with Gaussian shape, normalized amplitude 0.5 and 
width £ rj 2 /im (so that x c = 0.5). The theoretical 
prediction is that such a soliton is trapped if its velocity 
is below 0.55c (because fi cr it ~ 0.14, see Figure|3J). How- 
ever, in this experiment, we can consider variations of the 
external magnetic field B around the value 352 G, where 
the scattering length has the minimal value rs —0.23 nm. 
Increasing the field to the value B = 450G we can in- 
crease the scattering length to the value w —0.18 nm. 
This means that the scattering length can be varied by 
25%, and thus a nonlinear barrier V n i with normalized 
amplitude 0.5 can be generated. The theoretical predic- 
tion is that, in the presence of the linear trap and the 
nonlinear barrier, the soliton will be trapped only if its 
velocity is below 0.2c (because fi cr u ~ 0.05, see Figure 
0}, and it will be transmitted otherwise. This means that 
the nonlinear modulation dramatically enhances the do- 
main of parameters for which solitons can be transmitted 
through the linear trap. 

V. CONCLUSION 

In this paper we have investigated the time-dependent 
nonlinear scattering of bright solitonic matter waves 
through different types of barriers. We have consid- 



ered the transmission of matter waves through inhomo- 
geneities in the form of localized linear and nonlinear 
potentials. The adiabatic dynamics of the wavepackets 
as well as the radiative processes during the transmission 
have been analyzed. 

To analyze the dynamics we use the perturbed 1ST the- 
ory, which allows us to predict the trapping of a bright 
soliton by a trap potential and the reflection of a soli- 
ton by a barrier potential. The parameters (mass and 
velocity) of the transmitted soliton can be estimated by 
computing the radiated mass and energy and by using 
the conservations of the total mass and energy. The en- 
hanced transmission of a soliton through a linear trap 
by a nonlinear modulation of the scattering length is ex- 
plained by this theory. 

The analytical predictions have been checked by com- 
parisons with numerical simulations of the GP equa- 
tion. The formulas are valid in the asymptotic frame- 
work where the potentials have small amplitudes. It 
turns out that, for small or moderate potential ampli- 
tudes |Vjo| < 0.5, \V„m\ < 0.5, the perturbed 1ST theory 
gives quantitatively accurate predictions. For large po- 
tential amplitudes |Vjo| > 1, \V n io\ > 1, the theoretical 
predictions of the perturbed 1ST theory are still in quali- 
tative agreement with the numerical results. This means 
that the perturbed 1ST theory is useful for probing the 
parameter space and exhibiting interesting phenomena. 
The enhanced transmission can be observed in the exper- 
iments with bright matter wave solitons in elongated trap 
with proper variation in space of the external magnetic 
field and the trap potential 0, 0] . One of the problems 
that should be addressed for future consideration by this 
approach is the nonlinear resonant scattering on the (pe- 
riodic or random) chain of nonlinear barriers fl9| and the 
scattering on time-dependent linear and nonlinear barri- 
ers [2(1 . The latter problem is important for many areas 
of condensed matter. 
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